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We generalize a positivity constraint derived initially for parity-conserving processes to the parity- 
(~| . violating ones, and use it to derive non-trivial bounds on several Sivers functions, entering in the 

CIh' theoretical description of single spin asymmetry for various processes. 
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^ ' The relevance of positivity in spin physics, which puts non-trivial model independent constraints 

. on spin observables has been already clearly demonstrated. Different methods can be used to establish 
these constraints and many interesting cases have been presented in a recent review article [Ij. As a 
Tij- ■ well known example let us remind the positivity bound on the quark transversity distribution 5q{x, Q"^) 
^T) ■ [1], which survives up to next-to-leading order QCD corrections. This distribution plays a crucial role 
I in transverse spin physics. In this letter we will show that one can also derive non-trivial positivity 
bounds on Sivers functions [3], an important physical quantity for the interpretation of transverse 
single spin asymmetries (SSAs). Let us first recall a powerful general result which constrains the SSA 
in an inclusive reaction of the type + p^ — )• C + X, where both initial protons are transversely 
^ I polarized and C is any particle produced at rapidity y, whose polarization is not analyzed. It has 
been proven [1] that for y = and parity-conserving processes 



i-^,vjv(y = 0) >2|Ajv(y = 0)| (1) 

where Ajsr^ is the double transverse spin asymmetry. In the following we will first generalize Eq. ([T]) 
to include parity-violating processes and show it is valid in full generality for both parity-conserving 
and violating ones, then use it in the case where C is to put non-trivial bounds on the Sivers 
functions. 

To start, we consider the inclusive reaction of the type 

A{spml/2) + B{spml/2)^C + X , (2) 

where both initial spin 1 /2 particles can be in any possible directions and no polarization is observed 
in the final state. The observables of this reaction, which are the spin-dependent differential cross 
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sections with respect to the momentum of C, can be expressed in terms of the discontinuities (with 
respect to the invariant mass squared of X) of the amphtudes for the forward three-body scattering 



A + B + C^A + B + C, (3) 

as given by the generahzed optical theorem. 

The spin-dependent cross section corresponding to ([2]) can be defined through the 4x4 cross 
section matrix M and the spin density matrix p as [H d] 

aiPa,Pb) = Tt{Mp), (4) 

where Pa, Pb are the spin unit vectors of A and B, p = pa ® Pb is the spin density matrix with 
Pa = {I2 + Pa ■ o^a)/2, and similar for p;,. Here I2 is the 2x2 unit matrix, and a = {crx,cry, Uz) stands 
for the 2x2 Pauli matrices. The cross section matrix M can be parametrized as 

M = do [h + AaNCTay ®l2 + Abj^h dby + A^NCTay <8) dby + AnCTaz CFbz 

ax ^ ^bz 

+AaLCrz (8) /2 + Ahlh <8) AaSCTx O /2 + Absh CTx 

+ALNcrz <^cry + A^LCTy <^ cr^ + A^sf^y <^ cr^ + Asnctx <^ CTy] , (5) 

where ctq stands for the spin-averaged cross section, and L, S, N are unit vectors, in the center-of-mass 
system, along the incident momentum, along the normal to the scattering plane which contains A, B 
and C, and along N x L, respectively. It is important to realize that for a parity-conserving process 
there are only eight independent spin-dependent observables [HllllS]: the unpolarized cross section, 
two single transverse spin asymmetries AaN and AbN, and five double spin asymmetries Ann, -^ll, 
-^ss, ^LS and AsL- Now for the most general case including also the parity-violating processes, one 
has sixteen independent spin-dependent observables. Besides those in the parity-conserving processes, 
one has four additional single spin asymmetries Aai, A^i, Aas and A^s-, and four additional double 
spin asymmetries Aln, Anl, A^s and Asn- 

Eq. ([5]) is fully justified, since substituting it back to Eq. (jH) one has 

a{Pa,Pb) = (yo[^ + AaNPay + AbNPby + ANNPayPby + ALLPazPbz + AssPaxPax 
+ALsPazPbx + AsiPaxPbz + AalPaz + AhlPbz + AasPax + AbgPbx 
+ALNPazPby + A^LPayPbz + ANsPayPbx + AgNPaxPby] (6) 

The crucial point is that M is a Hermitian and positive matrix and in order to derive the positivity 
conditions, one should write the explicit expression of M as given by Eq. ([5]). In the transverse basis 
where ay is diagonal, we have found that the diagonal matrix elements Ma are given by 

Mil = {1 + Ann) + (AaN + AbN), (7) 

M22 = (1 - Ann) + {AaN - AbN) , (8) 

M33 = {I - Ann) - {AaN - AbN) , (9) 

M44 = {1 + Ann)- {AaN + AbN). (10) 

Since one of the necessary conditions for a Hermitian matrix to be positive definite is that all the 
diagonal matrix elements has to be positive Ma > 0, we thus derive 

l±ANN>\AaN±AbN\. (11) 
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Since both Ai\f and ^atat are parity-conserving asymmetries, it is not surprizing that this positivity 
bound obtained for parity-conserving processes is preserved. 

Back to the case for + — )• C + X where the initial particles are identical, we have AaN{y) = 
—AijN{—y)- Using this relation in Eq. ([IT]) , one obtains for y = 0, 

l-ANNiy = 0)>2\ANiy = 0)\. (12) 

We thus found Eq. ([1]) is valid in full generality, for both parity-conserving and parity-violating pro- 
cesses. In fact we found that physically 1 — Ajstn i: 2Ai\f at y = is nothing but the cross sections 
C7(ti)/o'o or a{\.'\')/ao at 90 degrees, which must be positive. Here t and \. denote transversity states. 

Now let us study the implication of Eq. ([T]) in a parity-violating process p^ + p^ — )• + X. 
The double transverse spin asymmetry A^jy for this process has been studied in Ref. [H [71 El [S]- 
The current conclusion is that AjVTV is expected to be negligible. The reasoning is the following: 
in collinear factorization formalism, the relevant function - the transversity distribution does not 
contribute as pointed out in [UlTj. This is due to the fact that the bosons are pure left-handed, 
their electroweak couplings do not allow a right-left interference. Higher-twist contributions have been 
studied in Ref. [8], where it was shown that a large suppression occurs from Sudakov factors, which 
increases with Q, the invariant mass. On the other hand, in the /cj^-factorization formalism, as shown 
in full details in Ref. [9] , one could (potentially) have additional non-zero TMD contributions El: a) 
contribution involving the product of hi and/or hyp, and b) contribution involving fyj, f^p and qit Qit- 
We want to point out the terms from (a) do not contribute for W production. This is because hi and 
h^P have the same Dirac structure as transversity in collinear factorization, thus the same reasoning 
for transversity not contributing in collinear factorization applies also here. On the other hand, for 
the terms from (b), there could be a Sudakov suppression \1Q\ applies to such terms that makes the 
contributions very small. Even if such a Sudakov suppression does not apply, according to current 
phenomenology since and qit are much smaller than /i, they are expected to generate very small 
effects to ^TVAT. Thus we will assume Ajqjy k, oH, and therefore Eq. ([1]) reduces to 

1>2|A^(2/ = 0)|, (13) 

to be compared with the usual trivial bound 1 > |^Ar(y = 0)|. There are also some studies that 
Aj^iN could receive non-zero contribution from the physics beyond the standard model [SKU], but we 
will not discuss this in our paper. We will now see how this strong result can be used to put severe 
constraints on the Sivers functions. 

Let us start with the exact definition of the quark Sivers function. Following the convention of 
Ref. |12] . the transverse momentum dependent (TMD) quark distribution in a transversely polarized 
hadron can be expanded as 

/,/;,t(x,k^,5) = f,,h{x,k^) + ^A^/g/^t(x, A:^) 5 • X kx) , (14) 

where p and k_|_ are the unit vectors of p and k_|_, respectively, fq/hi^^ ^±) is the spin-averaged TMD 
distribution, and A'^ fg/f^T-{x, k±) is the Sivers function [3j. There is a trivial positivity bound for the 
Sivers functions which reads 

\A^fg/f,,{x,k^)\<2fg/h{x,k±). (15) 

^Note in Ref. [9], these terms are written for Drell-Yan dilepton production, not for W production. 
•^If ^jvjv is found to be somewhat much larger than our expectation in the future, our positivity bounds could be 
modified accordingly using Eq. (|12|l . 

'*For non-trivial bounds involving several different TMD distributions see [16] and also Eq. (4.86) in Ref. [1] 
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Let's see how the SSA of the W bosons could give further non-trivial bounds on the Sivers functions. 
The SSA of the W bosons in hadronic collisions A^B W + X has been studied in terms of the TMD 
factorization formalism in [13\ 114], It is given by the ratio of the spin-dependent and spin-averaged 
cross sections, 



A 



da 



N 



dAa{S±) 



(16) 



where y and are the rapidity and the transverse momentum of the W boson, respectively. The 
spin-dependent and spin-averaged cross sections are given by 



dAa{S±) 
dyd^cij_ 

da 
dy (i2q_L 



2 



a,b 



S± • {pa X K±)^^ fa/At{Xa,ka±)fb/B{xb,kb±) 
CrO \^ab\'^ J fa/A{Xa, ka±)fb/Bixb, h±) , 



(17) 
(18) 



a,b 



where the simple integral symbol represents a complicated integral, / = J d'^]i.a±d'^'kb±6'^ {q^± — k^x — 
kb_|_). The partonic cross section uo is given by ctq = (7r/3)\/2 GpAf^/s with the Fermi weak coupling 
constant Gp and s = {pA +Pb)'^- J2ab runs over all light (anti)quark flavors. Vat are the CKM matrix 
elements for the weak interaction. Xa and Xb are the parton momentum fractions given by 



Xn 



Mw 

7^ 



Xb 



Mw 
7^ 



(19) 



To further simplify the expression of Ajy, it is usually assumed in the phenomenological studies 
that the x and dependence of the TMD distributions can be further factorized as follows \12\ [T5] 



fq/hix,k. 



A^fg/,,,{x,k^) 



Mx)gik^) , (20) 
A''f^/^,{x)h{kMk±), (21) 
For the k± dependence, a Gaussian ansatz is usually introduced [El [15], which has the form 

1 



g{k^) 



Akl) 

Y^k^^-ki/M! 

Ml 



(22) 
(23) 



The smart choice of h{k±) is introduced in [12], which satisfies h{kj_) < 1. Thus the positivity bound 
in Eq. (jlSp implies 



A^/,/,t(x) <2/,(x). 



(24) 



With the Gaussian ansatz, one could carry out the integration in Eqs. (jlTl) and (jTSl) analytically and 
obtains 



AN{y = 0) 



Eab\Vab\'^''fa/Ar{x)fb{x) 



Y.ablVabl'^ faix) Mx) 

where x = M-^ j ^/s for y = from Eq. ([TUj) . and H(q±) is given by 



(25) 



Hiq±) = S± ■ (pa X q_L) 



2e 



2\2 



Ml [{kl) + (fc2)]2 



<fcl)-(fc|> 



(26) 
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where {k^) = {k'j_)/[Mi+{k'j_)]. The g_L-dependent function H{q±) reaches its maximum H{qj 
when ql = {kl){{kl) + {k^,)/{{kl) - {k'^^)), with H{q^)^^, given by 



H{qA 



2\2 



[{kl) + {kIW 



2{kl) jkl) + (fe,^) 
{kl)-{kl) 



(27) 



Using the fact that 1 > 2|j4iv(?/ = 0)| for any q± and ^/s in Eq. ([25]) . we thus derive a new bound for 
the Sivers functions 



J2ab\yab\'^ fa{x) fbix) 

For W~^, Eq. (p8]) can be simphfied as 



< 



1/2 



H{q±)r 



(28) 



u(x) 



(i(x) 



< 



H{q 



max 



(29) 



While for W , one obtains the following constraint 



d{x) 



u{x) 



< 



H{qi_\ 



(30) 



Eqs. (j28p - (j30p are the main results of our paper Even though we derive the bounds for the Sivers 
functions in production in pp collisions (or Drell-Yan (DY) type process), our results also apply 
to the Sivers functions in semi-inclusive deep inelastic scattering (SIDIS) process. This is because one 
could show from the parity and time-reversal invariance that the Sivers function in SIDIS and in DY 
have the same functional form but opposite signs |17] . 

If we further use the values for the parameters and (/c^) from the global fitting of the experi- 
mental data for the SSA of the hadron in SIDIS processes, say, from Ref. |12j . 



{k] 



0.25 GeV^ Ml = 0.34 GeV^ 



(31) 



is M? 



one obtains H{qi_)^a_^ k, 0.31, thus 1/-H^((7±)max ~ 3.2. The fitted value of in Ref. [12j is ivi^ 
0.34to;?6. If one uses = 0.64, we get a better bound, 2.6 instead of 3.2. And at the same time the 
Sivers functions get bigger as shown in Ref. [12j. The constraints become better. While the previous 
bound in Eq. p5p would give the number 4 in the right hand side of the inequalities in Eqs. ()29p and 
(|30p . thus the new bound is slightly stronger than the previous one. However if one uses the lower 
value Ml = 0.18, the bound becomes 4.75, which is larger than the trivial one. Also if one increases 
[kj) while keeping the same M^, the bound will become looser. For example, for M^ = 0.64, when 



2 \ > 0.75 GeV^ 



then the bound will become looser than the trivial one. 



{k 

To summarize, in the framework of the factorized Gaussian ansatz, we have derived non-trivial 
bounds for the Sivers functions associated with n, d, u, and d quarks. The bound is slightly stronger 
than the previously derived one on the Sivers functions. Our results impose important constraints 
in the global analysis of the Sivers functions in the phenomenological studies of the single transverse 
spin asymmetry. We emphasize that the significance of our results should not be restricted to one 
special parametrization based on a very limited set of data. It is also important to realize that our 
newly-derived bounds Eqs. (f28]) - p0]) can be used as a consistent check of the quark Sivers functions 
extracted from the global analysis of the single transverse spin asymmetries. 



^There are some sign errors in Ref. [13], which have been corrected in the l.h.s of Eqs. (|29l [30 
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Finally this also applies to the gluon Sivers function which is accessible in the SSA for direct photon 
production [18|, because in this case also it was found that A^ttv is very small |19| [20], say hardly 
more than 1 or 2%. The most appropriate reaction is in fact pp^ — )• 7 + jet + X, which will deserve 
further considerations. We look forward to future accurate measurements of single spin asymmetries 
to greatly improve our knowledge about Sivers functions. 

J.S. acknowledges some interesting discussions with X. Artru. We thank A. Metz for helpful discussions 
and useful comments. Z.K. is grateful to RIKEN, Brookhaven National Laboratory, and the U.S. 
Department of Energy (Contract No. DE-AC02-98CH10886) for supporting this work. 
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